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This  study  was  carried  out  in  an  effort  to  determine  the 
longitudinal  stability  derivatives  of  a  Ryan  Mavion  airplane  using  dynamic 
flight  testing  methods.  This  was  done  by  recording  the  transient  response 
of  the  airplane  to  elevator  forcing  functions  in  the  form  of  steps  or 
pulses.   Data  reduction  was  carried  out  using  the  equations  of  motion 
method  and  the  derivative  method.   Due  to  partial  failure  of  the  instru- 
mentation the  results  from  the  equations  of  motion  method  were  disappointing. 
More  successful  results  were  obtained  using  the  derivative  method  in  con- 
junction with  the  useful  data  acquired  from  the  equation  of  motion  method. 

Concurrent  with  this  report  an  investigation  was  carried  out  to 
find  the  stability  derivatives  of  the  same  airplane  from  both  theoretical 
considerations  and  from  steady  state  flight  testing.  This  investigation 
is  described  in  Reference  7.   An  effort  has  been  made  to  correlate  the 
flight  test  results  obtained  from  the  dynamic  and  steady  state  methods  and 
predictions  made  from  theory.  This  correlation  appears  in  Appendix  B. 

Reasonable  values  of  C^   were  obtained  from  dynamic  testing 
which  were  consistent  with  speed,  suid  position  of  center  of  gravity  changes. 

Dynamic  testing  also  gave  reasonable  values  for  the  damping 
derivatives  which  were  consistent  at  the  sane  speed,  but  which  indicated 
an  increase  in  the  damping  derivatives  with  increase  in  speed. 

The  proper  changes  in  values  for  the  stick  fixed  msineuver  margin 
were  found  as  the  position  of  the  center  of  gravity  was  changed,  however 
a  decrease  in  maneuver  margin  was  found  as  speed  was  increased. 

The  elevator  power  was  found  with  good  consistency  as  the  center 
of  gravity  was  changed,  but  decreased  in  value  as  the  speed  was  increased. 
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All  values  of   Cm  r  were  somewhat  high  as  compared  with  theoretical  values. 

An  attempt  was  made  to  explain  tho  changes  in  Cm r  and  stick 
fixed  maneuver  margin  that  occurred  with  speed  changes  in  terms  of  aircraft 
non-linearities. 
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INTRODUCTION 

rVnamic  testing  of  aircraft  for  stability  parameters  has  received 
increased  emphasis  with  the  advent  of  guided  missiles,  sonic,  and  super- 
sonic aircraft.  Some  of  the  major  efforts  in  this  field  have  been  dynamic 
testing  of  a  B-25  by  Marshall  E.  Mulllns  at  Edwards  Air  Force  Base,  (Ref. 
4) I  the  dynamic  testing  of  an  F-30A  by  R.  C.  Kidder  at  Cornell  Aeronautical 
Laboratories,  Inc.  (Ref.  5)i  and  guided  missile  testing  by  The  Naval  Air 
Missile  Test  Center,  Pt.  Mugu,  California.   In  all  previous  cases  no  effort 
has  been  made  to  conpare  the  results  with  those  obtained  by  steady  state 
flight  test  methods  for  the  same  aircraft.  This  report  is  part  of  a 
triple  project  on  a  Ryan  Navion  in  which  dynamic  results  are  to  be  compared 
with  theoretical  and  steady  state  results.  The  actual  comparisons  are 
included  in  Appendix  B. 

The  dynamic  testing  with  which  this  report  is  concerned  consists 
of  analysis  of  the  transient  response  of  the  aircraft  to  an  elevator  forcing 
function.  Using  various  forms  of  the  equations  of  motion  the  stability 
parameters  for  the  airplane  were  reduced  from  the  data  by  the  method  of 
averages,  least  squares,  and  by  the  use  of  graphical  plots. 
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THEOHI  AND  ANALYSIS 

The  normal  procedure  in  the  determination  of  aircraft  transient 
motions  is  to  solve  the  differential  equations  which  describe  the  aircraft 
motion.   In  doing  this  the  constant  coefficients  in  the  equations  of  motion 
are  estimated  from  theoretical  considerations  or  measured  by   steady  state 
test  methods. 

The  problem  dealt  with  in  this  report  is  the  inverse  of  that 
stated  abovei  ie.,  the  determination  of  the  coefficients  in  the  equations 
of  motion  if  the  time  history  of  the  aircraft  transient  motions  are  known. 
The  response  history  of  the  aircraft  can   oe  displayed  as  a  function  of 
frequency  as  well  as  time.  In  theory  the  response  in  either  the  time  or 
frequency  danain  can  be  analyzed  to  determine  the  constant  coefficients 
in  the  equations  of  motion.  This  report  deals  with  the  determination  of 
the  constant  coefficients  of  the  longitudinal  equations  of  motion  from 
transient  responses  displayed  in  the  time  domain. 

Development  of  Equations: 

Throughout  this  report  it  is  assumed  that  the  incremental  changes 
in  forwsird  velocity  are  negligible.  With  this  assumption  the  equations 
of  motion  as  developed  in  Ref.  1  aret 


^(X  ^  olc<  -ci©  =  0  (1) 


c.o(^c..do(^  Cr.,je^i^c^^e^cS--o 


d5^^  <^    ^  (2) 
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If  elevator  lift  is  consif'ered  Important,  its  effect  can  be  inclu'ied  in 
the  lift  equation  since; 

^  -  C  mj- 

"-h       '  ^  (3) 

Substitution  in  equation  (1)  yields: 

Solving  equation  (4)  for  olo(  and  substitution  in  equation  (2)  yields » 


(5) 


Noma!  acceleration,  r» ,  and  elevator  angle,  6   ,  may  be  substituted  for  (X 
since: 

Substitution  of  this  expression  for  0(  into  equation  (5)  yields: 


<%-^Ty    ^  (^-.e^^-Jde^^^.K^^J^-^d'^^O    (^) 
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This  equation  if  integrated  once  with  respect  to  y^.    yields! 


~  V\(^e  -Q ^  (3) 


\ 


It  is  possible  also  to  eliminate  all  variables  except  a),  6     and  their 
derivatives  from  equations  (2)  and  (4).  Solving  equation  (4)  for  dO  and 
d  6      and  substituting  these  into  equation  (2)  yields: 

•^^      ^^/cy^        ^^/c  (9) 

Solving  equation  (6)  for  do(  ,  and  do(  in  terns  of  n  and  ^   yields: 

^^-      Cl,%  (U) 

Substituting  equations  (6),  (10),  and  (11)  into  equation  (9)  yields: 
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^^cX       ,    O^  \         C  /  C, 


^^^..  ^^^;/7  =^ 


Equation  (12)  if  integrated  once  may  be  expressed  as: 


(12) 


^"(tf- '"¥)!"  -  ''-^(t-  H)P '  -      ■"> 


where 


v^  f   ^^  -  ^-</v-^-<^. 


>)n-  ^  dS  ^C^(c.,^.C.^:jS-^  h  J 


n 
(14) 


Throughout  the  remainder  of  this  report  the  following  quantities 
are  defined  as: 


A  = 
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C  mn  ^  (~  rTt  cj  e( 


-Qi-^-  ^^  (1=) 


■o( 


C^^^    ^      Cr^j^  (16) 


C^J/-/-        ^^^     1  (17) 

Ct,  (  £::^  ~f  Qr^)  (i3) 


M   = 


Substituting  equation  (lo)   through  (19)   into  equations   (7),    (3)   anri   (13) 
yields: 


A  n    -i-  dcIO  ^  C  S  ■=   hJ  ^  (20) 


^n     -^  SS      ^  CfS    -   /^  J  6 


n    •+• 


'S=  V 


(21) 
(22) 


Equation  (22)  may  be  further  rearranged  as: 

.        ^  ^  -  ^  ^  -  / 

Solution  of  Equations: 

In  order  to  solve  the  equations  of  motion  for  values  of  the  deriva- 
tives, one  quantity  must  be  known.   In  this  case  the  known  quantity  was  the 
airplanes  moment  of  inertia  about  the  Y  axis.  To  determine  the  moment  of 
inertia,  the  airplane  was  oscillated  about  the  jack  points  as  shown  on  Fig.  1. 
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If  it  is  assumed  that  the  airplane  as  suspended  represents  a  second  order 
system  with  zero  damping,  it  is  possible  to  deterraine  the  moment  of  inertia 
about  the  jack  points,  from  the  formula  shown  on  Fig.  1,  The  moment  of 
inertia  was  then  transferred  to  the  center  of  gravity.  The  position  of 
the  center  of  gravity  was  determined  by  using  scales  under  each  jack  point 
and  under  the  tail  for  two  different  inclinations  of  the  airplane.  With 
the  moment  of  inertia  about  the  center  of  gravity  known  it  was  possible  to 
determine  /'   about  the  center  of  gravity.  By  keeping  an  accurate  record 
of  weight  added  to  and  removed  from  the  airplane  ^l    was  determined  for  each 
run.  The  accurate  cietemination  of  this  quantity  for  each  run  is  very 
important  as  small  errors  in  A  cause  large  errors  in  the  solution  for  the 
constants  in  the  equations  of  motion,  equations  (20),  (21),  or  (23). 

Examination  of  the  equations  of  motion  shows  it  would  be  most 

desirable  to  measure  <3e  ,  '^  »  c/(9  ,  and  a  O    directly.   Since  an  angular 

2 

accelerometer  was  not  available,  direct  measurement  of  d  9  was  impossible. 

Equation  (6)  shows  the  desirability  for  having  normal  acceleration  data  in 

order  to  determine  Ci^   .  For  this  reason,  and  realizing  that  there  are 

certain  inherent  difficulties  associated  with  the  accurate  determination  of 

a,  it  was  decided  to  use  nonnal  acceleration  as  one  of  the  measured 

variables.   In  order  to  make  the  same  instrumentation  system  available  for 

use  with  steady  state  data  and  so  that  the  elevator  hinge  moment  equation 

might  be  considered  in  future  studies,  it  was  decided  to  include  stick 

force  as  one  of  the  measured  variables.  The  resulting  instrumentation 

yielded  simultaneous  values  of  a,  n,  ^p    ,  dG,  and  F  •  The  actual  circuits 

^  s 

used  are  described  in  Appendix  A. 
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Since  h  is  a  measured  quantity,  measurement  of  n,  S^  ,  d0,  and 

2 

d  9  at  three  different  values  of  time  would  permit  the  solution  of  equation 

(20). for  A,  B,  and  C.   In  like  manner  measurement  of  j  n,  6, /'^^  aiid  d0 
at  three  values  of  time  would  permit  the  solution  of  equation  (21)  for 
A,  B»  and  C.  Knowledge  of  j  n   and  j  ^^    ,  and  v  for  two  values  of  time 
would  permit  the  solution  of  equation  (22)  for  L  and  M« 


Methods  of  handling  redundant  data: 

If  the  data  were  perfect,  evaluation  of  the  variables  at  only- 
three  insteuits  of  tine  would  be  sufficient  to  uniquely  determine  the 
coefficients  of  equations  (20)  and  (21).  With  less  than  perfect  data, 
however,  some  method  must  be  used  to  determine  the  set  of  coefficients  which 
most  logically  fit  the  data.  The  authors  have  used  several  methods  of 
handling  the  redundancy  arising  from  evaluation  of  the  variables  at  many 
values  of  time.  The  first  was  a  simple  application  of  least  squares  assuming 
the  error  to  be  all  in  one  variable!  the  second  was  a  method  of  averages: 
the  third  was  a  graphical  plotting  method.  This  latter  method  is  the  most 
desirable  since  it  gives  a  physical  picture  of  the  data  scatter  and  the 
fairing  process.  The  lack  of  physical  feeling  for  the  fairing  process 
when  least  squares  or  averages  are  used  to  reduce  the  data  is  of  great 
importance . 

When  the  method  of  least  squares  was  used  the  normal  equations  of 

2 

equation  (20)  assuming  all  errors  in  d  6  are: 

A^n"-  +  8^-» da  +  C^nZ  '  A  s:-?c <y fe  (24) 
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In  applying  the  method  of  least  squares  It  was  necessary  to  carry  as  many 
places  as  possible  In  ccsnbining  the  variables  and  solving  the  simultaneous 
equations.   No  justification  is  offered  for  the  apparent  inconsistency 
of  carrying  eight  to  ten  significant  figures  in  the  summations  of  the 
various  combinations  of  variables  when  the  variables  themselves  cannot  be 
taken  from  data  sheets  with  more  than  three  significant  figures.   It  is 
merely  stated  as  a  finding  that  as  many  places  as  possible  must  be  used 
in  order  to  get  reasonable  answers  using  least  squares  procedures.   It 
follows  that  least  squares  procedures  should  not  be  attempted  without  the 
use  of  an  eight  t6  ten  place  calculating  machine.   It  is  also  to  be  noted 
as  a  corollary  of  the  above  that  least  squares  procedures  are  very  unfor- 
giving of  small  errors  in  data  transcription,  summations,  etc. 

The  method  of  averages  consists  of  simply  using  the  average  values 
of  n,  d0,  S   ,  and  d  9  taken  from  random  groupings  of  the  individual  data. 
The  best  results  wei^  obtained  by  simply  grouping  the  data  into  approximately 
even  thirds.  The  first  of  the  three  simultaneous  equations  then  would  be 
made  up  of  the  summations  of  the  variables  over  the  first  third  of  the  time 
traces  considered  appropriate.  This  method  consisted  of  dividing  the  data 
into  three  such  groups.  The  simplicity  and  speed  are  such  great  advantages 
over  even  the  simplest  application  of  least  squares  that  its  applicability 
should  be  at  least  evaluated  in  any  extensive  data  reduction  problems. 

Practical. Considerations 

Data  for  eight  runs  were  reduced  using  any  or  several  of  the 
methods  described  above.  The  results  were  consistent  for  different  methods 
applied  to  the  same  run  but  the  only  answer  which  seemed  appropriate  was 
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the  value  for  B  =  ^^^a  ^^An^'     "^^^  values  of  A  emd  C  were  of  improper 
magnitude  and  inconsistent  from  run  to  run.  V/ith  this  in  mind  the  equations 
for  the  unknowns  A,  B,  C  in  detenninant  form  were  examined. 


/I     = 


9 

Zd9 

z.j 

/) 

td'-^ 

ZJ6 

ZJ 

fc^'a 

Ide 

2S 

In 

Fc/« 

^S\ 

In 

'^J9 

ZcT 

In 

Ide 

^J 

(2  7) 
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In 
In 

^  n- 


1  Ad. 


IS 


IS 


(8S) 
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In 
In 
In 


I6e        Id'-O 
Ids       Zd  'a 

Icie     76'e 


In  Ide  ^^ 
In  Ic^B  IS 
^n        ^^e       ^S 
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Exaniination  of  equations  (27),  (23)»  and  (29)  shows  that  equation 
(23)  is  the  only  one  where  an   error  in  the  pitch  rate  gyro  calibration 
would  not  affect  the  results. 

An  examination  of  the  pitch  rate  gyro  revealed  that  considerable 
friction  had  developed  in  the  bearings  and  it  was  very  doubtful  if  the 
gyro  wheel  was  getting  up  to  governing  speed. 

If  it  is  assumed  that  the  gyro  maintains  a  constant  speed  during 
the  1.5  seconds  during  which  data  is  being  taken  for  any  particular  run, 
the  reduced  values  of  B  will  be  correct.  Values  of  A  and  C  will  of  course 
be  directly  dependent  on  the  gyro  speed  for  that  particular  run. 

Because  of  tiir.e  limitations  it  was  impossible  to  repair  the  gyro, 
recalibrate,  and  take  new  data.  This  necessitated  re-examination  of  the 
equations  of  motion  in  order  to  eliminate  pitch  rate  as  an  unknown.  The 
method  is  very  similar  to  that  proposed  in  Ref.  2  and  called  the  derivative 
method. 

An  examination  of  equation  (14)  shows  that  if  (C//?^^  +  C„,^^   )  and 
C^^  are  known,  an  estimation  of  C y^j ^     will  define  the  quantity  zf    . 
Equation  (23)  is  then  recognized  as  that  of  a  straight  line  where  the 
variables  are  ■%4/'  and  ■'ylr  and  the  intercepts  with  the  coordinate  axes 
determine  values  of  L  and  M.   At  each  instant  of  time  it  is  only  necessary 
to  compute  /^^  ,  f^     ,   and  V   i  make  a  plot  of  equation  (13)  on  linear 
coordinates!  and  fair  the  best  straight  line  through  the  points  to  determine 
the  intercepts.  These  graphical  plots  were  extremely  enlightening.  They 
indicated  that  L  can  be  determined  to  much  greater  accuracy  than  can  M. 
A  blind  application  of  least  squares  or  averages  could  never  indicate  this 
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conc^ition.  Knowledge  of  this  condition  would  be  valuable  for  deciding  n^ 
which  variable  to  miniitize  the  error  for  application  of  least  squares. 

It  should  also  be  noted  that  the  answers  are  very  insensitive  to 
the  estimated  values  used  since  the  known  value  of  {C/y,,  +  C/^,  ,^)  makes  up 
the  major  part  of  v. 

In  the  process  of  data  reduction  it  was  attempted,  with  poor 

results,  to  use  a  data  instead  of  n  for  solving  the  moment  equation.  The 

fifteen  hundred  series  nans  yielded  very  poor  results.  This  was  accounted 

for  by  the  fact  that  an  accurate  in-flight  calibration  of  the  angle  of 

attack  indicator  to  account  for  vane  position  error  effects  was  never 

attained  for  these  runs.  The  nature  of  the  difficulty  is  described  in 

Appendix  A. 

The  a  data  was  used  in  equation  (6)  to  obtain  values  of  C.  .  The 

a 
use  of  a  data  necessitates  taking  into  account  the  effect  of  up  wash  and 

the  effect  of  pitching  velocity.  The  effect  of  up  wash  is  determined  from 
steady  flight  calibration  of  the  angle  of  attack  indicator.  The  effect 
of  pitching  velocity  is  determined  by  measuring  the  distance  of  the  vane 
ahead  of  the  center  of  gravity  and  computing  the  velocity  ccmponent  at  the 
vane  due  to  pitch  rate.  Another  effect  that  it  was  not  possible  to  account 
for  quantitatively  was  the  effect  of  acceleration  in  pitch.  Due  to  flex- 
ibility of  the  boon  angular  acceleration  of  the  airplane  caused  bending  of 
the  boom  in  such  a  meinner  that  individual  values  of  the  C^   determined  at 
each  instant  of  tirae  yielded  a  plot  as  shown  below.  Better  stiffening  of 
the  boom  should  eliminate  this  source  of  error. 

As  the  dynamic  response  of  the  angle  of  attack  vane  was  not  deter- 
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rained,  any  natural  dynamic  response  of  the  vane  consisted  of  an   error  in 

the  a  data  of  undetermined  magnitude. 

For  the  above  reasons  a  data  was  used  only  to  detenr.ine  C.   from 

a 
the  sixteen  hundred  runs. 


Q. 


Zf/M( 


The  solution  of  the  equations  of  motion  by  any  of  the  techniques 
used  does  not  allow  solution  for  unique  values  of  the  derivatives  ^^7i>:joc  » 
Cmac      *     ^^'S     *  ^"^  ^^(i&    '     ^^  order  that  the  results  may  contain 
approximate  values  for  the  above  derivatives  it  was  assumed  that  C'/n^^  J^^^j^ 
and  —^      =   .5.   It  is  important  to  realize  that  without  the  above  assumption 

Ofoc 

dynamic  flight  test  procedures  as  described  herein  do  not  yield  unique 
answers  for  the  individual  derivatives. 
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PROCEDURE 

The  procedure  followed  in  this  report  may  be  considered  in  the 
following  parts:  1.  Instrumentation,  2.  Calibration,  3.  Flight  test, 
4.  Data  reduction. 

1.  Instrumentation 

The  design,  construction,  and  installation  of  the  instruments 
used  in  this  report  took  approximately  four  months.  The  circuits  actually 
used  as  well  as  the  design  criteria  may  be  found  in  Appendix  I.  The 
reasons  for  the  choice  of  variables  to  be  measured  may  be  found  in  the 
Theory  and  Analysis  section  of  this  report. 

2.  Calibration 

The  calibration  curves  for  each  circuit  used  may  also  be  found  in 
Appendix  A.   A  brief  explanation  of  inherent  difficulties  in  measuring 
certain  variables  may  be  found  in  the  Theoi^  and  Analysis  section  of  this 
report.   Calibration  was,  of  course,  carried  out  continuously  as  the  flight 
test  progressed  to  insure  that  the  calibration  curves  remained  applicable. 
The  one  exception  to  the  above  statement  was  the  calibration  of  the  pitch 
rate  gyro.   It  was  calibrated  in  the  final  stages  of  instrument  installation 
at  the  Naval  Air  Test  Center,  Patuxent  River,  Maryland,  and  no  further 
calibrations  were  performed.  The  difficulties  encountered  in  the  early 
stages  of  data  reduction  emphasized  the  desirability  of  continuous  calibra- 
tion checking. 

3.  Flight  test 

The  flight  procedure  for  transient  dynamic  testing  is  extremely 


-15- 

siraple  amd  represents  the  greatest  advantage  of  dynamic  testing  over 
steady  state  testing  nethods.   Smooth  air  is  essential  in  order  that  a 
good  steady  state  flight  condition  may  be  attained.   After  attainment  of 
a  good  steady  state  condition,  the  steady  state  is  disturbed  in  some  pre- 
determined manner.  Pulse  type  elevator  forcing  functions  were  applied  by 
simply  moving  the  yoke  a  predetermined  amount,  holding  it  in  the  new  posi- 
tion for  approximately  one  second,  and  then  returning  it  to  the  original 
position.  The  resulting  pulses  were  poor  approximations  to  an   analytical 
pulse.  However,  in  an  equations  of  motion  technique  this  was  not  important. 
The  best  method  found  for  applying  a  step  elevator  input  consisted  of 
holding  the  stick  against  the  force  caused  by  a  preset  value  of  trim  tab 
and  then  siiriply  letting  go  of  the  stick.  The  elevator  is  so  heavily 
damped  that  the  resulting  elevator  motion  was  very  nearly  a  step. 

All  data  was  taken  at  19.5  inches  manifold  pressure  and  about  1850 
engine  revolutions  per  minute.   Low  engine  speeds  were  required  to  reduce 
angle  of  attack  boom  resonance.  Variations  in  speed  were  obtained  by 
varying  the  rate  of  descent.   Data  was  taken  for  two  center  of  gravity 
positions  and  two  speeds  using  both  pulse  and  step  elevator  forcing  functions. 
4.  Data  reduction 

The  various  methods  of  data  reduction  are  described  in  the  Theory 
and  Analysis  section  of  this  report. 
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SAMPLE  CALCULATIONS 

A  sample  of  each  type  of  calculation  is  carried  out  in  this  section 
to  further  clarify  the  procedure  used.  The  calculations  are  carried  out 
on  run  fl'1640.  The  values  of  the  variables  taken  from  run  #1640  may  be 
found  in  Table  II.  Table  I  contains  the  various  constants  required  in  the 
calculations  and  the  calibration  factors  taken  from  the  calibration  curves. 

Equation  (2l)   is 

A/n    -^  S  e   ^    C/J   =    ^  do 

The  values  for  Jn   «  6  ,  /cT  anj  J^may  be  obtained  for  each  point  indicated 
in  Table  II  since: 

where  C>^  -  calibration  factor  for  n  in  c|'s//^. 

JSdPr   =     yI'^'^^    ^'^'^^    seconj.) 
Jded'^/l-    =    ^    ^     Qq    fedi     ^^nc^  seconds) 


de  = 


do 


d'/r   '    '^^    '     T  Q,^   r/ncAes) 


With  values  of  the  variables  at  each  point  three  simultaneous 
equations  may  be  obtained  as  follows  where  the  summations  are  over  the 
various  points  of  timet 
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'^%Sn%    *BZ&    ^     Cij^d\=    hide 
Aiifodi'^    ^3i&     ^    CZf'id*''7   =    Aide 


Actual  calculations  yielded  the  following  simultaneous  equations: 

-  .000091  A  +  .017371  B  -  .020173  C  =  .010376 

-  .198996  A  +  .110355  B  -  .074817  C  =  .020030 

-  .601369  A  +  .204940  B  -  .125027  C  =  .014376 

Solution  of  these  equations  simultaneously  for  B  gives  B  =  -  .141341. 

It  is  clear  that  equation  (7)  as  it  stands  could  be  solved  using 
values  of  J  ,  n,  d6,  and  d  9.  These  variables  in  inches  and  seconds  units 
are  also  included  in  Table  II. 

If  equation  (6)  is  rearranged  it  is  possible  to  solve  for  C,  . 

'  r    -        -  Cm  ^  .^/  ^ 

o( 

Values  of  n,  ^  ^   and  o<  may  be  obtained  from  Table  II  using  the  calibration 
factors  of  Table  I.  The  value  of  0(   so  obtained  should  be  further  cor- 
rected for  the  error  introduced  in  the  vane  indicator  by  pitching  velocity. 

r 


A  o(    ^       ^ 


1/ 
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where  7  ia  true  velocity  in  ft/sec. 

It  is  necessary  to  estimate  C^  ^   .  ^//ic   "^s  estimate^  theoretically  to 

be  -1.3.   An  incorrect  value  of  C^^.  ^  vill  of  course  be  reflected  in  an 

incorrect  value  of  C,  but  the  percentage  error  in  C   will  be  much  less 

a  a 

than  the  percentage  error  in  C  ^^  ^  .  Using  all  points  of  Table  II  the 

following  equation  was  obtained  : 

i2 


S37 


L  r^      - 


S3  7 


--  S.OC 


For  this  particular  case  an  increase  in  ^  y,i  c     of  bO^>   gives  an 

increase  of  12$   in  C.  . 

a 

The  terras  on  the  right  hand  side  of  equation  (13)  have  been  esti- 
mated as  shown  at  the  end  of  this  section  for  runs  #1522,  1523,  1640,  and 
1652.  These  estimated  quantities  have  been  determined  using  values  of 
{Crri      "'"  ^'>^d&^   determined  by  equations  of  motion  method  and  as  rationalized 
in  Fig.  2.  An  average  value  of  C^ ^   and  the  theoretical  value  of  C-ni$~  "/^S 
have  also  been  used  in  these  estimations.  Tables  III,  IV,  V,  and  VI  contain 
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the  calculations  for  v  anr*  the  variables  "^y^X  and  /(f   .  These  variables 
are  plotted  and  the  values  of  L  and  M  determined  on  Figs.  3  and  4. 
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RESULTS 

The  values  of  C.      were  consistent  with  theoretical  considerations 
and  no  trends  were  indicated  as  speed  and  center  of  gravity  position  were 
changed. 

The  damping  derivatives  were  found  with  good  consistency  at  the 
sarae  speed  but  an  increase  in  damping  derivatives  (  C^^   +  6^.  )  was  found 
with  increase  in  speed. 

r  Cr,     r      1 

The  values  of  stick  fixed  maneuver  margin  ^  -—3   /   ^rC&j 

indicated  a  change  of  the  proper  amount  when  center  of  gravity  position 

was  changed  at  constant  speed.  However,  a  decrease  in  stick  fixed  maneuver 

margin  was  found  as  speed  was  increased.  The  stick  fixed  maneuver  point 

for  90  mph  was  58^,  for  110  mph  it  was  525^. 

The  value  of  (_^        was  found  with  some  consistency  at  each  speed, 

o 
but  a  decrease  was  indicated  in  C,^     as  soeed  was  increased.  The  values 

of   C^     were  somewhat  high  in  comparison  with  theoretical  values. 

o 

The  only  immediately  apparent  enigmas  in  the  results  then  seemed 
to  arise  as  a  result  of  speed  changes.  This  would  lead  one  to  suspect  an 
error  in  some  or  all  of  the  constants  used  in  data  reduction  which  were 
functions  of  speed.  These  constants  are  ^  ,  and  C  .  The  values  of  2^   and 
C  used  were  very  carefully  examined  and  no  errors  could  be  found. 

Li 

During  the  in-flight  calibration  of  the  angle  of  attack  indicator 
a  condition  was  found  which  might  explain  some  of  these  changes  with  speed. 
A  sharp  increase  in  up  wash  at  the  wing  tip  was  noted  as  speed  was  decreased 
below  35  raph  indicated.  The  test  aircraft  was  equipped  with  a  spoiler  near 
the  wing  root.  The  possibility  existed,  then,  that  non-linearities  in  flow 
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a.nr'   lift  were  occurring  during  pull  ups  froT  the  90  mph  steady  state. 

If  non-linearities  in  the  flow  di'^  exist,  they  would  be  reflected 
in  low  values  of  normal  acceleration.  Further  the  actual  error  in  the 
normal  acceleration  would  become  greater  with  time.  The  recorded  values 
of  nonnal  acceleration  would  be  less  by  a  progressively  larger  amount  than 
that  which  would  have  been  recorded  sans  non-linearities.  The  recorded 
values  oT     \ /?   would  have  been  less  but  not  by  as  large  a  percentage  as 
normal  acceleration  itself.  The  error  in  aC/f    was  not  as  apparent  as  that 
for  n  and  j  n   . 

The  equation  from  which   cL   an"^  stick  fixed  maneuver  margin  were 

s 

obtained  was: 

Z    %      .    ^  ^    =   / 

The  quantity  v  as  determined  in  Tables  III,  IV,  V,  and  VI  was  predominantly 

made  up  of  n,  particularly  after  the  first  few  tenths  of  a  second.  The 

1/1 
non-linearities  discussed  then,  if  present,  would  have  tended  to  make   ^ — 

and  ^  ^     too  large,  the  error  increasing  with  time.   Since  ^ was  a 

larger  quantity  than   -i^   these  errors  would  tend  to  decrease  the  slope 

of  the  lines  faired  into  the  data  in  Figs.  3  and  4.  This  could  have  given 

the  increased  apparent  stability  found  at  lower  speeds.  Also  since  the 

points  at  later  times  were  used  predominantly  in  fairing  Figs.  3  and  4, 

the  increased  value  of  ^      at  90  mph  could  slLso  be  caused  by  these  non- 

linearities. 

An  attempt  to  evaluate  the  effects  of  non-linearities  on  the  values 

of  the  damping  derivatives  as  deterriiined  by  equations  of  motion  method   lead 
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to  no  (definite  conclusions.   It  is  possible  that  some  of  the  apparent 
increase  in  damping  derivatives  with  speed  may  have  been  caused  by  low 
values  of  normal  acceleration  at  the  lower  speed.  The  increase  in  damping 
derivatives  with  speed  indicated  an  increase  in  tail  efficiency.  The  test 
procedure  required  as  low  an  R^A   as  possible  to  reduce  the  angle  of  attack 
boom  resonance.   Power  used  was  19.5  Inches,  full  low  RRn  at  both  speeds. 
This  required  an  increase  in  rate  of  descent  with  increase  in  speed.  This 
does  not  appear  consistent  with  the  increase  in  tail  efficiency  unless  the 
tail  was  operating  in  a  particularly  unfortunate  position  with  respect  to 
wing  wake  at  the  lower  speed. 
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C0NCLDSI0N5  AND  RECCMTffiNDAT  IONS 

"tynamic  flight  testing  is  in  general  a  perfectly  feasible  method 
of  accurately  determining  the  stability  derivatives  of  an  airplane. 
Certain  basic  requirements  for  obtaining  good  results  are  as  follows: 

1.  Inclusion  of  pitch  rate  as  one  of  the  measured  variables. 

2.  Weight,  moment  of  inertia,  and  airspeed  must  be  known  very 
accurately. 

3.  Instrumentation  must  be  very  accurate  without  unknown  instru- 
ment dynamic  characteristics. 

Dynamic  testing  has  the  following  advantages  and  disadvantages. 
Advantages: 

1.  Simplicity  of  flight  procedure. 

2.  Short  duration  of  in  flight  testing. 

3.  Non-linear  effects  may  be  included  in  the  equations  of  motion. 

4.  Values  of  the  derivatives  obtained  are  the  ones  the  airplane 
feels  in  response  to  transient  maneuvers.  They  are  the  values  for  velocity 
equal  a  constant. 

5.  Close  correlation  exists  between  frequency  response  deterrr.ined 
from  dynamic  testing  and  auto  pilot  design  requirements. 
Disadvantages: 

1.  Instrumentation  is  generally  very   expensive. 

2.  Data  reduction  is  lonf^,  tedious,  and  subject  to  error. 

3.  Some  of  the  results  can  not  be  compared  to  static  flight  test 

2 
results  which  are  for  C^V  equal  a  constant. 

In  recommending  that  the  work  be  continued  the  most  important 
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improvements  in  the  methods  used  are  listed  below. 

1.  Further  efforts  should  be  made  to  obtain  accurate  pitch  rate 
data. 

2.  A  pendulum  or  turntable  should  be  constructed  for  calibration 

of  the  rate  gyro  and  accelerometer. 

2 

3.  An  angular  accelerometer  would  permit  the  measurement  of  d  9 

directly. 

4.  Angle  of  attach  data  would  be  improved  by  boom  stiffening, 
and  dynamic  calibration  of  the  vane. 

5.  A  frequency  analysis  of  the  airplane  response  might  yield 
interesting  results. 
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APPENDK  A 
Description  of  Test  Equipment 

The  test  airplane  was  a  standard  Ryan  Navion  instrumented  to  record 
time  histories  of  elevator  angle,  angle  of  attack,  stick  force,  pitch  rate 
and  normal  acceleration.  These  variables  were  recorded  with  a  Consolidated 
Engineering  Corporation  type  5-116-P3-14  recording  oscillograph.  A  wiring 
diagram  for  the  instrumentation  system  \i3e6   is  shown  in  Fig.  5.   Schematic 
diagrams  showing  the  operation  of  the  individual  circuits  appear  as  Figs. 
6  and  7. 

Certain  basic  considerations  were  comon  to  the  design  and  operation 
of  all  of  the  circuits.  These  were: 

(a)  To  make  the  equivalent  resistance  of  each  circuit  be  350  ohms 
as  seen  from  the  oscillograph  leads,  in  order  to  meet  the  manufacturer's 
recanmendation  for  62^   of  critical  damping  for  optimum  galvanometer  response. 

(b)  To  operate  the  sensitivity  adjustment  of  each  circuit  at 
about  50fc  attenuation. 

(c)  To  use  as  small  balancing  resistors  as  possible  in  order  to 
maintain  linearity  of  the  output,  and  at  the  same  time  keep  the  resistors 
large  enough  to  limit  the  current  drawn  from  the  battery. 

(d)  To  obtain  maximum  galvanometer  deflection  for  the  range  of 
variables  measured  at  each  flight  condition.  The  galvanometer  deflection 
was  one  inch  per  12.4  M  amps  of  current. 

Power  for  the  instruments  was  obtained  from  a  45  volt  dry  cell  battery 
with  a  center  tap.  Voltage  from  this  battery  to  the  bus  lines  was  regulated 
to  20  volts  for  all  tests  using  a  variable  resistor  and  voltmeter. 
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Slevator  An^le  and  Anple  of  Attack  Circuits.  The  ^    and  of  circuits 
are  shown  schematically  in  Fig,  7.  The  potentiometers  used  to  sense  the 
amgular  motions  were  12,500  ohn,  Micro-}^ax  precision  potentiometers  with 
outputs  guaranteed  to  be  linear  within  one  per  cent.   The  motion  of  the  cc 
vane  and  the  elevator  was  transmitted  to  the  potentiometer  shafts  in  a  one 
to  one  ratio.  The  2500  ohm  balancing  potentiometer  is  used  to  adjust  the 
circuit  output  to  zero  at  various  elevator  or  a-vane  positionG.  The  addi- 
tional 5900  ohn  resistors  in  series  with  the  balancing  potentiometer  were 
required  to  limit  the  current  drawn  from  the  dry  cell.  The  1500  ohm  load 
resistor  was  required  to  adjust  the  current  to  the  galvanometer  to  the 
desired  magnitude.   Sensitivity  adjustment  was  achieved  through  a  500  ohm 
impedance  pad.  The  1000  ohm  variable  shunt  resistor  was  placed  in  parallel 
with  the  circuit  to  make  the  galvanometer  external  damping  resistance  be 
350  ohms  as  required. 

The  calibration  of  the  elevator  angle  circuit  is  shown-  in  Figs.  8 
and  9.   These  calibrations  were  made  using  a  propeller  protractor  mounted 
on  the  elevator  to  measure  A  ^^j  .  Fine  adjustments  of  the  elevator 
position  were  made  by  moving  the  elevator  with  a  hydraulic  jack  under  the 
trailing  edge. 

Calibration  of  the  angle  of  attack  circuit  was  made  in  flight  using 
the  propeller  protractor  to  measure   A  "^  j"  ,  Calibrations  were  also  made 
on  the  ground  using  a  specially  made  fitting  on  the  boom  to  measure  the 
angle  between  the  vane  position  and  the  boom  centerline.  These  calibrations 
are  shown  in  Figs.  10  and  11,  suid  Illustrate  the  difference  in  calibration 
under  flight  conditions  caused  by  airloads  on  the  vane  and  boom.  Due  to 


\ 
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an  oscillograph  breakdown  it  was  impossible  to  get  an  in-flight  calibration 

for  the  1500  series  runs  before  the  sensitivity  and  zero  balance  settings 

were  changed.  Therefore  the  in-flight  calibration  for  the  1500  series  runs 

was  made  using  the  difference  between  ground  and   flight  calibrations  for 

the  1600  series  runs  and  extrapolating  for  a  correction  to  be  added  to  the 

1500  series  ground  calibration. 

Stick  Force,  Pitch  Rate  and  Normal  Acceleration  Circuits.  The  oper- 

ation  of  the  F  ,  9.  and  n  circuits  are  similar  and  are  shown  schematicsdly 
s  •' 

in  Fig.  6.  The  strain  gauge  bridge  circuits  used  to  measure  these  variables 
will  be  described  individually  later.  The  5000  or  50,000  ohm  balancing 
potentiometer  in  each  circuit  gave  a  zero  output  adjustment  to  the  circuit 

« 

over  a  wide  range  of  values  of  F  ,  G,  and  n.  The  5000  or  50,000  ohm  re- 
sistors  in  series  with  the  balancing  potentiometers  were  necessary  to  orovide 
a  current  liir.it  to  protect  the  galvanometer  in  case  the  circuit  should 
inadvertently  be  turned  on  with  the  sensitivity  adjustment  full  open  and 
the  balancing  potentiometer  turned  against  the  stops  in  either  direction. 
The  sensitivity  adjustment  and  damping  resistors  are  the  same  as  described 
for  the  a  and  ^^   circuits. 

The  accelerometer  used  in  the  n  circuit  was  a  model  F-1. 5-330  accel- 
erometer  manufactured  by  the  Stathara  Laboratories,  Incorporated.  The  manu- 
facturer' s  specifications  required  an  13  volt  input  to  the  acceleroneter 
which  was  obtained  with  a  100  ohm  variable  resistor  in  series  with  the 
circuit  leads  to  the  power  supply.  The  accelerometer  had  a  range  of  +1.5  g'  s. 
The  circuit  was  calibrated  between  zero  and  one  g  sjid  assumed  to  be  linear 
to  1.5  g.  The  calibration  curve  is  shown  in  Fig.  12.   The  noise  level  in 
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the  n  circuit  was  extremely  sensitive  to  airplane  vibrations,  particularly 
to  engine  rpm.   In  order  to  minimize  this  the  7r-section  filter  shown  in 
Fig.  6  was  U3e(^  to  filter  out  vibrations  above  12  cycles  per  second. 

The  stick  force  bridge  was  made  up  of  3R-4  strain  gauges  cemented 
to  the  spokes  of  a  specially  constructed  force  wheel  which  replaced  the 
standard  pilots  wheel  in  the  airplane.  Calibration  of  the  stick  force 
circuit  is  shown  in  Fig.  13. 

The  pitch  rate  circuit  bridge  was  made  up  of  3R-4  strain  gauges 
mounted  on  a  leaf  spring  which  restrained  the  precession  of  a  gyroscope 
when  the  gyroscope  was  subjected  to  an  angular  velocity.  The  3R-4  strain 
gauge  bridge  requires  a  15  volt  input  which  was  obtained  by  using  a  100  ohm 
variable  resistor  in  series  with  the  circuit  leads  to  the  power  supply. 
The  pitch  rate  gyro  was  calibrated  on  a  calibrating  pendulum  at  the  Naval 
Air  Test  Center,  Patuxent  River,  Md.,  and  subsequent  recalibrations  were 
not  feasible.  This  calibration  established  the  linearity  of  the  circuit, 
and  is  shown  in  Fig.  14.  This  linearity  extended  throiagh  six  points  not 
appearing  within  the  range  of  Fig.  14.  A  250,000  ohm  ±1.1%   and  a  50,000  ohm 
±•1.1%   precision  resistor  was  placed  in  parallel  with  one  side  of  the  bridge 
as  shown  in  Figs.  5  and  6.  At  the  time  of  the  original  pendulum  calibra- 
tion the  bridge  unbalance  caused  by  these  calibration  resistors  was  noted 
to  give  a  circuit  output  equivalent  to  particular  values  of  pitch  rate. 
Then  with  the  instrument  installed  in  the  airplane  a  calibration  could  be 
made  for  each  test  flight.  This  was  done  by  momentarily  closing  a  calibra- 
tion resistor  switch  with  the  aircraft  in  steady  flight.  A  galvanometer 
deflection  was  thus  obtained  for  a  known  value  of  pitch  rate.  This  point 
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along  with  the  point  of  zero  deflection  for  zero  pitch  rate  was  then  used 
to  establish  a  calibration  curve  for  the  flight. 

The  electric  motor  in  the  gyroscope  was  powereH  from  a  24  volt 
storage  battery.   It  was  discovered  during  the  data  reduction  that  the 

# 

gyro  motor  had  had  a  progressive  bearing  failure,  and  that  0  values  were 
unreliable.  This  failure  could  have  been  discovered  had  a  calibrating 
pendulum  been  available  for  calibration  after  each  test  flight. 
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APPENDIX  A 


The  following  Is  a  tabulation  of  the  results  from 
Princeton  Reports  No.  231  and  232.  Report  No.  232  consists 
of  the  detennination  of  longitudinal  stability  parameters, 
for  a  standard  Ryan  Navion,  using  theoretical  calculations 
and  steady  state  flight  test  techniques.  Report  No.  231 
consists  of  determining  as  many  as  possible  of  these  same 
parameto? s  for  the  same  airplane  using  dynamic  flight  test- 
ing techniques. 
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Clean  -  Power  Off 




Parame  ter 
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.090 

=m. 
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ZxpjnXn'^.fAon  of  ths   t?,ble3   ".bovo   nhoiJG   the   folio:, 'Ins: 
1.   '^aluoc   of  C/^  aotemined  bv   t-ho   differont  nethodc 
'.-.'orG  all  -.rit'.iin  12,t, 

?•      Theoretic".!   r.nd   itc-  ly  stato  v.lnoa   of  C7ng''.z^''^Q 

v\-^r-j  clorol:''  '.;?]-. le   tho  v.luo  of^^^fror.  dyn-^/^io   to^tin^  zee 
to  bo   ?bout  301  hi3h.      Somo  of  this  coMld  be   iccoLmted  for 
b^^  tho   fact  th-^.t  no-'.i-lino'^.rities     in  the   olov"tor  ^,n3le 
ooter.tionotor  u-iod  for   .Ivnainio  ---^^"z  o?.U3e>'    -     ^  —  -'i '^i  "> -'.t;'^ 


o 


f  10-15"'  error  in  Idividual  vv.luoc   of  th^levator  riii^lO' 


/5 

'■^o'joTor,    it^lnorob.?.ble   that  the   Blopo   of  tha  elevator  --^.n^^le 
c"..li' ••'•  ■^:* -^  ^   c^irvo   could  be  noro   thin  10"'  off.      In  pXcd'j 
st'^te   toiitinc*    oo-jer  for  level  fli'jht  v.'!'.3  used.      In  the 
djnarp.ic   testing,   lor>s   po-.-er  ••r?.s  nned  in  ordor  to  ]ceep  the 
c-i:ine  R-"-  lo^r  to   -void      -n^le   of   -tticl:  boor-    r'-'^on.'^.nce. 
The   110  ;  /H  rir  noeod  vras  att?.ined  by  eatVolishin^/  a  ote-idy 
rate   of  descent. 

3,    "^'^-^   valiiec   of  the   ntic':  fi::od  iTTorMn/'^r     ■'^int     /^»7y^/c 
a^ain  '\j;reo   bot-.^een  theoretical  ind  steady   otate  valueo  vrhile 
t'-^e  value  fron  dynnjnic   tsntins  ^eenc   to  be   17>^  hi^h.      The 
dlfforenoe   reems  exconr.ive   since   It   re  :;  re  cents   a  47"^  deocrep- 
ancy   in  the   maneuver  r-^.r^^in.      The    graphical    olots  of  FI^g. 
3   and  ^  of  Pa  .ort  "o.    ^^31   indicate   the  diffic^ilty  of     obtain- 
ing froTi  dynamic  data  tho   intercept  v^hich  dotor''in^   the  "lan- 
ouver  mar::in.       ^ven  thou~l''   the  Intercept  is  defined  by  lines 
crosain:^  at   -^    call   an.~le   this   c-^n  not  •account  for     the    47;1 
di3orep?ncy.      Dynamic   •  or'c  v\s   severely  '^----•rod  by  failure 
of  the   r-ite   T/ro   nj\d  continued  inveatir. 'tion  -.dth     accurate 
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r^.te  3yro  ner.ouremonta  ma:'  J'leld  bstt&r  rer/.ilts.  Tho  effect 
of  G?.rryln3  the  lo-.:er  oo^rer   for  the  dyngjilc  te3tln3  is  alco 
app?.rent  vfhen  tho  po/or-off  mp.neuver  uaroln  froa  oter.d:/'  stnte 
resulta  is  considered.   If  the  values  f2>Dni  dyn-inlc  testing 
are  conpared  with  the  pov/er-off  v?.lu3s  of  ate-^.dy  sto.te  results 
tho  di::2repfinoy  is  ?7'^  indicating  thr.t  the  effect  of  carryin3 
thD  lo-'er  jo-.:er  in  dyngjiic  testing  may.  be  very  Ir^rje, 

The  v.lue  of  Cw^  oht^lnod.  fron  stc^.dy  st-^.te  te:;tin3 
is  ruoh  lo  .-or  th-^n  th«t  fron  dyn?inic  t^atin^  for  the  sarae 
rer.sons   causing  the  discrop^,ncy  in  the  'rianeuvor  mp.r^in. 
Tho  diffv^rence  in   <-?^^  is  further  ?,ccentu?.ted  by  the 
sniall  value  of  the  dampin.'^  deriv--tlvec  obt?,ined  from  ste-T'.dy 
state  testin3.   The  theoreticnl  v^lue  of   6w^  dotemined 
froEi  the  stick  fixed  noutr'^1  point  agrees  vrith  the  v^^.lue 
of  (^^   from  ste-jdy  state  tenting  p:nd  tho  theoretical 
value  of   (^^  determined  from  stick  fi"ed  nianeuver  point 
a3r9es  faily  v;ell  vith  the  v:^.lue  of  Cm^   determined  frnorr, 
dynamic  tests.   This  sho^:s  thr.t  dynanic  test  results  are 
very  closely  related  to  maneuvering  st-r.bility  vrhile  the 
stead:/  r:ta*e  results  are  irore  nearly  re  1  tod  to  ':tatic 
st-bility.   • 

The  values  of  the  danping  derivatives  ^>fr^^    '^"^"-  ^j^ 
abtained  fron  theory  and  d3mr:Tiio  testing  a^roe  closely  Thile 
those  obtained  from  steady  state  t  '^tin^  are  t^>o  lov%   Tliis 
can  be  orimarily  accounted  for  by  the  f'^ct  th"t  the  d-^jno- 
in3  derivatives  are  very  dlffic'ilt  to  o-tract  fron  the 
steady  state  test  data  ?-nd  are  oreaumed  to-  be  in  error. 
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4.      The   .'^rontest  vcrl."ncG  'bot-^een  t-ho  thoorotlcpl  and 
otead^'-  stite   fll::ht  to  at  result  s,   with  the  e-coeptlon     of 
par"',r.eterr;     .rovl-^unl^r  dlacuscod,   -."^.3   found  to  be   in  the 
Gtlcl:  froG  uicnouvQV  oolnt  for  the   po-er-on,   cloan  condi- 
tion,  r.nd  in  the  ctlck  fixed  neutr-^.l  .:^oint  for  the  pover- 
off,l^.ndi-'  condition   . 

It  is  b.lievad  th^.t   the  diacreoancy  in  the   oticl:  free 
maneuver   ooint  a'^.y  be   attributed  to  the  f-'.ct  th'^t  '^,t  lovr 
nor"-^!   '"ccol^^r••',tion3   the   stick  forco   re-^uAred  i-^.   lo-f  -^nd 
erroneous   indl c^.tions   of  ctick  force   are   poscible  bec?.uio 
of  the   friction  in  the   control  column  and  elevator  system, 
'Jn.tp.  rodM?tion   rnd  ?.nn.lviis   indlc^.ted  no   chin'je  in  nanou" 
vsrin^     oint     bet-.^^oen  nomal   n.c c e lor '^.t ions  of  1.0  to  1.3» 
■^ovrever,    if   ^n  intGrool?.tion  is  mr/Ie  betv:oen  the   ??tick 
free  rr'^.neuvor  -point   for  nom-^.l   ''.c colorations   of  1»0  to 
1.5,    ^.  viiue   of  'r-!,neuvor   point   of   ?._pro:;im'^.tely  5r\  •  me-n 
a9rod3''naniic   chord  is  obt-iined  uhich  corresponds- qnite   closely 
-.rith  the   53  '  nie?,n  .•■'.erodyn''jnic  chord  obt'^.lned   from   theoret- 
icr.l  conEid3r'^.tions» 

The   stick  fixed  neutr'?.!  point   for  the  pov;er-off ,   I'^.nd- 
in3     oonlition  v-hlch  V7.?,s   obt?.lned  from   f^teady   st'^te   to,:ts 
corres..ond3  to   the   trends   indlCT'.tQd  in  ti^e  other  po-.:er  con- 
oitions   rjid  confi-urations.      Ho;rever,   the  theoretic-.l  ina- 
lysis   iniic':'ted  consi^lor^bly  nor^   stability.      As'^il-rnnent 
of  t'^.c  error  involved  can  not  T?e  O.ono  vrith  '^,ny   Ao^ree  of 
'iccur'^.cy;  'ho'-rever,    it  Is  bellev.'^d  the  m.rjor  portion  of  the 
dlfferons-   co.n  be   '^ccoM.r.ted  for  in   tho   ovlu^tion  of/.f ".ntor 
(/-^^  -  ^Ajl  1  in  the   thoorotiG?.l   malysis. 
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